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Atmospheric Neutrino Flux: A Review of Calculations 

T.K. Gaisser^* 

'^Bartol Research Institute, University of Delaware, 
Newark, DE 19716, USA 

Interpretation and understanding of the evidence for neutrino oscillations depends on knowledge of the atmo- 
spheric neutrino beam. In this talk I assess how well various features are known. The goal is to determine to 
what extent uncertainties in the neutrino beam may limit the conclusions about neutrino properties and which 
features of the evidence for neutrino oscillations are most robust. 



1. Introduction 

The concept of using the atmospheric neutrino 
beam to look for neutrino oscillations is illus- 
trated in Fig. 1 §. With a single detector it is 
possible simultaneously to cover a range of path- 
lengths from ~ 10 to ~ 10^ km, corresponding re- 
spectively to downward moving and upward mov- 
ing neutrinos. The atmospheric neutrino beam 
has an energy spectrum determined by the steeply 
falling primary cosmic-ray spectrum, which gen- 
erates the neutrinos by interactions of the cosmic 
ray nucleons in the atmosphere. Examples of cal- 
culated neutrino fluxes are shown in Fig. 2. 

The neutrino spectrum, which falls with en- 
ergy, must be folded with the rising neutrino cross 
section to obtain the expected event rate as a 
function of neutrino energy. Most neutrino in- 
teractions are in the range from a few hundred 
MeV to a few tens of GeV. Thus atmospheric 
neutrinos have the potential to probe the range 
1 < L/E < 3 X 10"* km/GcV. From the standard 
two-flavor oscillation equation. 
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we find, therefore, that the atmospheric neutrino 
beam can in principle probe down to 6m? as small 
as ~ 2 X 10~^ eV^. (The actual lower hmit to 
sensitivity in 5rn? is somewhat higher than this 
because of the relatively large scattering angle be- 



tween neutrino and lepton in charged current in- 
teractions of low energy neutrinos.) 

Atmospheric neutrinos originate with the tt 
fj, —f e decay chain. Therefore at sufficiently low 
energy such that muons as well as pions decay, 
one expects 
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The anomalous value of this ratio [^|-|| (for which 
many sources of uncertainty cancel in the cal- 
culations) suggests neutrino oscillations involv- 
ing i-'p and/or v^. as a possible explanation. The 
telltale evidence Q for oscillations as the source 
of this anomaly comes from the pathlength de- 
pendence of the neutrino fluxes. The energy- 
dependence of the up-down asymmetry for muon- 
like events, coupled with the fact that events ini- 
tiated by electron neutrinos appear to have the 
expected energy and angular dependence, indi- 
cates that the primary effect is oscillations in- 
volving muon neutrinos with large mixing and 
« 3 X 10-3 eV^ g. 
In the remainder of this paper I will describe 
how the geomagnetic field affects the interpre- 
tation of the atmospheric neutrino data, review 
the primary cosmic-ray spectrum, discuss the un- 
certainties in the treatment of pion production 
and mention the comparison of the calculations to 
measurements of muons high in the atmosphere. 
It is important to review these points at this time 
because of the recent publication of two new cal- 
culations 0,^, which are three-dimensional, as 
compared to previous one-dimensional calcula- 
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Figure 1. Illustration of up-down symmetry of 
atmospheric neutrinos in the absence of effects 
of the geomagnetic field and assuming no oscilla- 
tions. 



tions. 

2. Geomagnetic field effects 

In the absence of oscillations, the only signifi- 
cant deviation from isotropy of the atmospheric 
neutrino beam is caused by the geomagnetic field, 
which prevents low energy primary cosmic rays 
from reaching the atmosphere at low geomagnetic 
latitudes. For example, at Kamioka the geomag- 
netic cutoff is sa 10 GeV for primary protons 
near the vertical. The downward neutrino flux 
at Kamioka is therefore lower than, for exam- 
ple, at Soudan, where the vertical cutoff is neg- 
ligible. Approximately half the downward neu- 



trinos with Ejy < I GeV at Soudan come from 
primary protons with i? < 10 GeV. This con- 
tribution is absent from the downward neutrino 
flux at Kamiokande. As a consequence the down- 
ward flux of sub- GeV neutrinos at Kamiokande 
is roughly half that at Soudan. The neutrino flux 
from the lower hemisphere is similar at the two 
detectors and intermediate between the down- 
ward fluxes at the two locations. This is because 
the neutrinos from below are produced over a 
large fraction of the Earth's surface and so av- 
erage over a range of high and low geomagnetic 
cutoffs. Thus the up-down asymmetry at each de- 
tector location is a combination of geomagnetic 
effects together with the effects of any oscilla- 
tions that may be present, and the combination 
depends on detector location. Moreover, the geo- 
magnetic effects become less important as energy 
increases. It is therefore useful to start by consid- 
ering the atmospheric neutrino analysis without 
the geomagnetic field. 

For a flux of primary cosmic rays that is spa- 
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tially isotropic, the atmosphere is a spherical sheU 
source of neutrinos with equal luminosity per unit 
volume independent of latitude and longitude. 
The number of neutrinos produced per unit vol- 
ume of atmosphere into solid angle d^l, 



s{e,h) 
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depends only on local zenith angle, 6, and alti- 
tude, h. In these circumstances one can show 
from the geometrical construction of Fig. 1, that 
the neutrino flux is up-down symmetric; i.e. sym- 
metric about cos 9 <-> —cosO. Since neutrinos 
from decay of pions, muons and other secondary 
cosmic rays are produced over a range of altitudes 
peaking around 15 to 20 km above sea level 
it follows that local variations in surface altitude 
introduce a negligible deviation from this symme- 
try. The symmetry also requires that differences 
caused by local variations of pressure are negligi- 
ble. 

A detector like Super-Kamiokande has an ac- 
ceptance that is up-down symmetric. In the ab- 
sence of geomagnetic effects, therefore, a sim- 
ple measurement of the up-down asymmetry as 
a function of energy is a probe of neutrinos oscil- 
lations. Therefore the simplest and most robust 
evidence for oscillations is a deviation from up- 
down symmetry in an energy range high enough 
so geomagnetic effects are small. Fig. 3 Q 
shows the distribution of neutrino energy for four 
classes of neutrino interactions. The correspond- 
ing distributions in energy per nucleon for the 
primary cosmic-rays are about a factor of ten 
higher. Thus the median primary of the multi- 
GeV events is about 50 GeV/nucleon, which is 
high enough so that geomagnetic effects are unim- 
portant. Neutrino-induced upward muons pro- 
vide a higher energy sample, but it is not pos- 
sible to make a simple up-down comparison be- 
cause the downward muon flux is dominated by 
muons produced from pion decay in the atmo- 
sphere rather than by interactions of neutrinos. 

Interpreting the up-down asymmetry of the 
neutrino fluxes in the GeV range and below re- 
quires that the geomagnetic effects be well under- 
stood. Geomagnetic effects on the primary cos- 
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Figure 3. Distribution of neutrino energies for 
various classes of atmospheric neutrino events. 



mic radiation are indeed very well understood. 
They form the basis for an entire subfield of 
cosmic-ray physics Low energy particles 

(few GeV) at low geomagnetic latitudes cannot 
reach the atmosphere to interact. Particles of in- 
termediate energy 10 to 20 GeV) that reach 
the atmosphere show a strong east-west asymme- 
try, while high energy particles 100 GeV and 
higher) are essentially unaffected by the geomag- 
netic field. 

The east-west effect is a consequence of bend- 
ing of positive primaries in the geomagnetic field. 
In fact, it was the observation of the excess of 
primary cosmic rays from the west p^Jl^ from 
which it could be deduced that the primaries were 
predominantly positively charged. The expected 
azimuthal dependence of the neutrino flux associ- 
ated with the east-west effect is the same whether 
or not the neutrinos oscillate because the distri- 
bution of pathlengths that contribute to a par- 
ticular zenith angle band is independent of az- 
imuth. Comparison between expected and ob- 
served azimuthal distributions is therefore a good 
check of the systematics of the whole chain of data 
analysis and calculations |l4] . The measurements 
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show the expected azimuthal dependence for both in Fig. 4. 

muon and electron neutrinos , indicating that 
the geomagnetic effects are weh understood. 



3. Primary spectrum 

A standard procedure for calculating the flux of 
atmospheric neutrinos is to generate atmospheric 
cascades for a spectrum of primary protons and 
nuclei and form the neutrino spectra in the ab- 
sence of the geomagnetic fields. The resulting 
neutrino spectra can be filtered through the ge- 
omagnetic configuration relevant for a particular 
detector location, discarding those neutrinos pro- 
duced by primaries that would not have reached 
the atmosphere. An alternative approach 
starts from the measured muon fluxes high in the 
atmosphere and uses the genetic relation between 
neutrinos and muons to obtain the muon spec- 
trum. In either case, the normalization of the 
neutrino flux depends on the absolute normaliza- 
tion of a measured spectrum of charged particles 
(either the primary cosmic rays or the muons). 
The advantage of starting with the muons is that 
one bypasses uncertainties in knowledge of pion 
production, which have a similar effect on both 
neutrinos and muons. A disadvantage is that 
the acceptance for muons is sensitive to details 
of propagation in the geomagnetic field. In ad- 
dition, the measurements of the primary flux are 
made with relatively long exposures at or above 
the atmosphere while the relevant muon measure- 
ments are generally made during a short balloon 
ascent. 

The mixture of nuclei in the primary spectrum 
is such that approximately 80% of nucleons in 
the cosmic radiation are free protons, 15% are 
bound in alpha particles and the remainder are 
in heavier nuclei. Spectra of protons and he- 
lium with i5 < 100 GeV/nucleon have been mea- 
sured at the top of the atmosphere in a series 
of balloon-borne spectrometer experiments. Re- 
cent measurements [^^-^ cluster around a lower 
normalization than an earlier standard refer- 
ence p3| . At higher energy, measurements have 
so far been possible only with calorimeters, which, 
because of punchthrough, may have larger sys- 
tematic uncertainties. The data are summarized 
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Figure 4. Primary spectra. 



4. Comparison of calculations 

There are now five independent calculations of 
the neutrino spectrum that start from the pri- 
mary cosmic-ray spectrum filtered by the geo- 
magnetic field. The calculations of Refs. 

are one-dimensional, assigning all produced 
neutrinos the directions of the primary particle 
that produced them. The ingredients of these 
calculations have been compared previously p^ . 
Recently two three-dimensional calculations have 
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been published . Treatment of the geomag- 
netic cutoffs, which now depend on 4 variables 
instead of 2, makes the calculation significantly 
more complex. These calculations are very im- 
portant because they check the major, technical 
simplifying assumption made by the previous cal- 
culations. 

A conclusion reported in Ref. is that dif- 
ferences for predicted event rates and how they 
depend on direction and energy are relatively 
small between one-dimensional and three dimen- 
sional versions of the same calculation. This has 
the important consequence that the simpler one- 
dimensional calculations can be used to explore 
the consequences of uncertainties in input to the 
calculations. Fig. 2 compares the Bartol neutrino 
flux with the 3-dimensional flux of Ref. 0. 
Both calculations use the primary spectrum of 
Ref. [Q, and the geomagnetic field has been 
turned off. Most of the difference, therefore, is 
presumed to be caused by differences in the treat- 
ment of pion production. 

5. Pion production 

Most pions with E < 100 GcV in the atmo- 
sphere decay before they interact. Therefore, for 
neutrinos in the sub- and multi-GeV range, only 
interactions of protons and helium play a sig- 
nificant role. The most important information 
needed about these interactions is the inclusive 
cross sections for pion production. The impor- 
tant range of interaction energies extends up to 
~ 100 GeV for sub-GeV and - 1 TeV for muhi- 
GeV interactions. The most probable energy of a 
primary proton for a sub-GeV event at Super-K is 
« 20 GeV. The corresponding number at Soudan 
is about 10 GeV because of the lower geomagnetic 
cutoff. The parent energies for multi-GeV events 
are correspondingly higher. The distributions of 
secondary nucleons is also important because a 
significant fraction of the neutrino production oc- 
curs in secondary or tertiary interactions of the 
nucleons. Neutral pions (and r] mesons) are also 
important in the sense that energy deposited in 
the electromagnetic part of the cascade is not 
available for production of neutrinos. Distribu- 
tions of kaons begin to be important for multi- 



GeV events (and they are dominant for neutrino- 
induced upward muons ||2^ ). 

Data on pion production have been discussed 
recently in Ref. [p^ . Existing measurements 
cover a significant fraction, but not all, of the 
relevant range of phase space for charged pion 
production in proton collisions on beryllium and 
aluminum. New, more precise measurements cov- 
ering all phase for proton interactions on a range 
of light nuclei (including nitrogen and oxygen) 
would be of great interest. Use of a helium beam 
would also be of interest. 

The difference between the neutrino calcula- 
tions shown in Fig. 2 most probably mainly re- 
flects a difference in the fraction of energy going 
into production of charged pions. The calcula- 
tion of Ref. j2^ uses a phenomenological hadronic 
event generator (TARGET) developed for cosmic- 
ray cascade calculations. FLUKA Q uses a 
more sophisticated microscopic model of particle 
production with intranuclear cascading. It incor- 
porates the event generator as an integral part 
of a cascade code capable of simulating interac- 
tions and cascades in complex detector geome- 
tries. The philosophies are quite different. The 
FLUKA interaction model is tested and adjusted 
by comparing directly to double differential cross 
sections for a wide range of data sets as measured 



{e.g. 



dpL d^pT 



The strategy with TARGET is to fit pr distribu- 
tions at each longitudinal momentum, p^, extrap- 
olate into unmeasured regions of phase space at 
each and integrate to obtain the energy flow 
into each secondary channel. 



dE^ 



A detailed investigation of the sources of differ- 
ence between these two models and others [p5| , 
and their implications, is currently in progress. 

6. Comparison to atmospheric muons 

Most muons and neutrinos are produced be- 
tween 10 and 30 kilometer altitudes |^ in closely 
related processes. Measurements of muons at 
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these altitudes therefore in principle provide a 
check of the neutrino calculations in which the un- 
certainties in pion production cancel to the extent 
that they are common to both the neutrinos and 
the muons. Recently there have been several com- 
prehensive measurements of muons during ascent 
of balloon payloads [ p^|p^j2^ , and some discrep- 
ancies with calculations have been noted. Gener- 
ally, the agreement is best at float altitude, where 
only the first interaction plays a role. An im- 
portant limitation to the use of muons to nor- 
malize the neutrinos, however, is that the muons 
are more sensitive to details of the calculation. 
For example, since most muons decay in flight, 
a change in interaction lengths, which moves the 
cascade up or down, can change the muon flux at 
a particular altitude without changing the cor- 
responding neutrino flux, which is an integral 
over the whole atmosphere. Since muons follow 
curved trajectories in the geomagnetic field, ap- 
proximating by straight lines can have a similar 
effect (by moving the muon decays lower in the 
atmosphere) . In addition, since positive and neg- 
ative muons have opposite curvature [pT| their re- 
sponse to cutoffs of the primary cosmic radiation 
are somewhat different. Three-dimensional calcu- 
lations of atmospheric muons are in progress. ]32| . 

7. Summary 

• The observed pathlength dependence of the 
muon-like events in Super-Kamiokande |^ 
points to neutrino oscillations as the source 
of the atmospheric neutrino anomaly. 

• Three-dimensional calculations of the atmo- 
spheric neutrino flux remove an important 
approximation present in previous calcula- 
tions. It appears, however, that the sim- 
pler one-dimensional calculations are ade- 
quate for exploring differences among cal- 
culations. 

• Uncertainty in the normalization of the pri- 
mary spectrum is now reduced to « ±15%, 
so the main remaining source of uncertainty 
is the representation of pion production. 

• If the neutrino flux is significantly lower 



than calculated in Refs. and |25[|, then 
there is a potential problem accounting for 
the relatively large number of electron-like 
events seen in SuperKamiokande within a 
predominantly ^ Vr oscillation scheme. 

• New measurements of atmospheric muons 
are being considered in airplanes and with 
slow balloon ascents. 

• Uncertainties in expected event rates due to 
the imprecise knowledge of neutrino cross 
sections in the GeV energy range have not 
been discussed here, but may be impor- 
tant, n 
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